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Onset of criticality and transport in a driven diffusive system 
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We study transport properties in a slowly driven diffusive system where the transport is externally 
controlled by a parameter p. Three types of behavior are found: For p < p' the system is not 
conducting at all. For intermediate p a finite fraction of the external excitations propagate through 
the system. Third, in the regime p > Pc the system becomes completely conducting. For all p > p' 
the system exhibits self-organized critical behavior. In the middle of this regime, at pc, the system 
undergoes a continuous phase transition described by critical exponents. 
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Cascade events where a single external perturbation 
may amplify in an avalanche like dynamics are seen in 
a variety of controllable systems, such as nuclear reac- 
tors O , avalanche diodes |0] and lasers [pj . In statistical 
physics one considers cascade or avalanche events in a 
certain type of slowly driven diffusive systems. The pro- 
totype of such dynamics is the sandpile model of Bak, 
Tang, and Wiesenfeld Q. This model exhibits large 
scale fluctuations with avalanches similar to the exam- 
ples mentioned above. The difference being that the 
slowly driven systems self organize to the critical transfer 
ratio, needed to see avalanches of all sizes; models ex- 
hibiting this phenomena are called self-organized critical 
(SOC). Following this idea, a number of similar models 
have been proposed. In particular, Zhang B and Manna 
1^ proposed sandpile models where the local redistribu- 
tion contained a random element. This defined another 
universality class than the original BTW model pj. Re- 
cently, similar ideas have been used in the explicit mod- 
elling (by MM ) of the observed rice pile dynamics in an 
experiment performed by the Oslo group [|lO|| . In fact, 
the rice pile models in MM belong to a universality class 
which encompasses a variety of different self-organized 
critical models, including both earthquake models, inter- 
face depinning models [Q and the Manna model p2[ . 
Discrete models in 1-d which do not fall into this class 
exist, namely the deterministic models in M] and their 
stochastic variants proposed in Q . 

In the present paper, we investigate the onset and ro- 
bustness of criticality in the above mentioned "rice-pile" 
class of models. To this end, we introduce a control pa- 
rameter p in a model which displays SOC, and investigate 
the robustness of the critical state as function of this con- 
trol parameter. The control parameter will influence the 
so called transfer ratio, J, which is defined as the proba- 
bility that an added particle will be transported through 
the system. When p is large the ratio will be critical, 
and when it is small the process will come to an end 
abruptly. The novel aspect of our study is that the model 
will be critical in a whole range of values of the control 
parameter. Further, in contrast to other SOC models 
with continuous parameters |l5[, the critical properties 



of our model does not depend on the control parameter 
as long as this is above a certain threshold point, p'. 

One may consider various formulations of stochastic 
SOC models in 1-dimension. One is the rice pile ver- 
sion where the particles perform a directed walk through 
the system, and another is the "slope version" where one 
only considers the equivalent redistribution of slopes (lo- 
cal stresses). As numerically demonstrated in 11^, the 
latter is similar to the Manna model Q . Here we present 
the model in the "rice-pile" language, since this allows 
for a direct counting of transport properties in terms of 
the number of particles which are absorbed/transported 
through the system. For physical applications one may 
consider the added particles as excitations. 

The model we study is defined as a simple extension 
of the one introduced in ||]. Consider a lattice [1, L] on 
which particles are added consecutively, to the column 
i — 1. The results to be reported in the following have 
been obtained for L — 200. The state of the system is 
defined in terms of the height profile h{i) of the pile, or 
equivalently in terms of the slopes Zi = h(i) — h(i + 1). 
The dynamics of the model is defined in two steps, of 
which 1) only takes place when 2) has come to an end: 
1) An avalanche is initiated by adding one particle to the 
column i = 1. If zi > 1 the column i = 1 is considered 
active. 2) During each step in the avalanche all columns i 
which are active can transfer one particle to the adjacent 
position i + 1 . The probability for each such toppling is p, 
and in case of a toppling at position z, then one particle 
is transferred from column i to i-\-l. For all toppled par- 
ticles, one tests whether any of the new local slopes zi-i, 
Zi, Zi+i fulfills Zj > 1. If this is the case for a column 
j, then one particle at column j is considered active in 
next step. The procedure is then repeated with the new 
active sites, until finally there are no more active sites. 
Then a new particle must be added to the first column 
{^ = 1). 

Notice that the model contains one parameter, p, 
which describes the toppling probability for slopes zi > 
Zc = 1. Thus, this toppling probability does not depend 
on the value of the slope when it exceeds Zc- As a result, 
the model allows for very steep slopes, and consequently 



also for very large storage of externally imposed particles. 
It is this open storage capacity which allows the model 
to display both subcritical and critical behavior. 

As usual, when considering systems that are slowly 
driven, one has to run the dynamics for a certain tran- 
sient time before its stationary features can be studied. 
In Fig. 1, we show the average transport ratio J, defined 
as the relative flux of outgoing particles to ingoing par- 
ticles in the stationary state. This is equivalent to the 
probability that an added particle at position i = 1 will 
escape at position i > L, as function of the parameter p. 
We observe three regimes: below p' « 0.53865 the system 
is nearly completely isolating and nothing penetrates. In 
the interval [p',pc], with pc « 0.7185, the system is par- 
tially conducting with a slope z = (zi) = [h{l) — h{L)]/L, 
which fluctuates around a constant value determined by 
p since both h{l) and h(L) grow at the same rate in this 
regime. Finally, for p > pc, the system is completely con- 
ducting, with a transport that exhibits the white noise 
behavior seen in Fig. 2. 

Before entering the detailed discussion of the behavior 
of the system around the two transition points, p' where 
conduction starts, and Pc where the motion of the pile 
gets pinned, we first investigate the fluctuations in the 
dynamics inside the pile. In Fig. 2 we show the fluctua- 
tions of 1 — J as function of time for two p values, p = 0.6 
and p = 0.75. For p < p' (and large L), no particles ever 
reach i = L (i.e., J = 0). For higher p values, the trans- 
port out of the system is noisy, with a behavior which 
changes dramatically when passing p = pc- In fact (see 
Figs. 2a and 2b) we would like to stress that although 
the dynamics in both intervals [p',Pc] and [pc, 1] are con- 
ducting, the fluctuations in the transport properties have 
a phase transition at Pc'. they exhibit a Brownian motion 
below Pc, but a white noise signal above. 

In Fig. 3 we plot the size of the relaxation events as 
measured by integrating over the number of topplings due 
to a single added particle. We see that for p > p' , the 
probability for having an avalanche of s topplings exhibit 
the power-law scaling 



p{s) oc l/s"'. 



1.57±0.05. 



(1) 



This result is in accordance with what we expect for the 
rice-pile universality class for 1-d systems driven at the 
boundary [||,|l|,|l|l. Thus, the SOC state in our sys- 
tem is not affected by the possibility that some particles 
permanently are absorbed in the system, as is the case 
for p' < p < Pc- The type of nonconservation imposed 
in this interval of p values does not destroy scaling. For 
p < p' the situation is different and the avalanche size 
distribution becomes exponentially bounded. When no 
particles are allowed to pass through the system, the sys- 
tem cannot build up long-range correlations. 

In Fig. 4, we show the system properties at respectively 
the onset to conductivity {p — p') and at the depinning 
transition (p = pc) ■ The onset to conductivity is the sim- 
plest to understand. For low p the chance that excitations 



propagate through the system is small. The cascade trig- 
gered by adding one particle to the pile at i = 1 dies out 
exponentially fast. Next, to understand the origin of the 
transition point p = p' consider the cascade process at a 
point i > 1. A particle toppling from position i makes 
particles at i — 1, i and i + l potentially active. If the cor- 
responding slopes take values larger than Zc{= 1), each of 
the columns will topple with probability p and make the 
new columns active. In particular, when part of the pile 
is in a state with slopes Zj ^ 1 then the toppling events 
will always generate new active columns which topple ex- 
actly with probability p. This means that the spreading 
of active particles are determined by directed percolation 
(DP) on a square lattice with 3 descendents. This process 
has a critical point at p = 0.53865 which is in complete 
agreement with our numerical estimate for p' , thus we 
use p' = 0.53865. Below the value p' the spreading of 
activity dies out and the Zi values increases indefinitely, 
making our subcritical DP picture self consistent. 

For p > p' , the DP becomes supercritical, and the 
system becomes at least partially conducting. Thus the 
local slopes do not grow indefinitely. Accordingly, the 
DP mapping breaks down, and the spreading of activity 
in the lattice will also be influenced by dynamically dis- 
tributed absorbing states (states with slope and 1) in 
the system. The transition ed p — p' may be quantified 
as in Fig. 4a, where we show the vanishing of the cur- 
rent as function of the distance to p'. In order to obtain 
a better scaling we subtract the finite current J' at p' , 
which occur only because we have a finite system. Given 
that p' = 0.53865 is exactly known from DP studies, and 
measuring the finite current for a system of size L — 200 
a,t p = p' to J' = J(0. 53865) = 0.013, we observe numer- 
ically that J{p) - J' oc (p - p'f with S' « 0.9. 

Approaching p = Pc, the system develops an ever in- 
creasing region extending from i = \ where many of the 
local slopes are small, hi p — pc, this region spans the 
entire system, and all particles added to i ~ 1 will even- 
tually be conducted out of the system at L (see Fig. 1). 
Above Pc the system profile remains bounded (pinned) at 
position i = L. By decreasing p, the profile starts moving 
upwards with a velocity that depends on pc — p. In Fig. 
4b, we see that the number of particles absorbed in the 
system per unit time scales as 



I ~ J cc\pc 



P\ 



6 = 0.9 ±0.1. 



(2) 



In the language of rice pile dynamics this would be the 
"velocity" of the pile (= h{L)/t), and accordingly the 
transition at p = Pc may be understood as a depinning 
transition for the profile of the rice pile. For p values on 
both side of this transition the system is critical, with 
exponents determined by the "Manna" universality class 
(with avalanche exponent r = 1.55 and avalanche dimen- 
sion D = 2.25 ± 0.05). It is an open question whether 
the depinning exponent 6 can be related to the universal 
exponents of this very robust class of SOC models. 



In summary, we have introduced and analysed a slowly 
driven diffusive system which exhibits two phase transi- 
tions: the first described by the critical point of directed 
percolation, and it defines a transition to a self-organized 
critical state described by a very robust universality class. 
The second transition, at pc, happens within the SOC 
state and can best be characterized as a transition to a 
state where the long-range meandering of the transport 
properties seen below Pc gets pinned at all values above 
Pc- We find it interesting that conductivity and the pres- 
ence of a SOC state is intimately related in the present 
modelling. 
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FIG. 1. Average transport ratio J of particles through the 
system, as function of the parameter p. Transport is defined 
in terms of fraction of conducted particles. Notice that when 
the average transport approaches unity, then in the rice-pile 
picture the growth velocity {— 1 — J) of the pile approaches 
zero. The point where the transport starts is p' ~ 0.53865, 
and the point where it reaches unity is pc ~ 0.7185. 
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FIG. 2. Temporal dependence of the average absorbtion 
rate (= "velocity" ) for (a) p — 0.6, and (b) p = 0.75. In 
both cases one observes fluctuations of about the same order 
of magnitudes. However, although the system is in the SOC 
state in both cases, the type of fluctuations are clearly differ- 
ent. For p > pc one observes Gaussian white noise, whereas 
p < Pc leads to long-range correlations (reminiscent of Brow- 
nian motion) in the transmitted signal. For p < p' everything 
gets absorbed, and there are no fluctuations. 
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FIG. 3. a) Log-log plot of the avalanche size distribu- 
tions for p — 0.4, p — 0.6, and p = 0.75. The observed 
power-law scalings correspond to the rice-pile universality 
class for p > p' , whereas the bounded avalanches for low p 
show that p < p' does not give critical fluctuations. 
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FIG. 4. Log-log plot of transport properties as a function 
of distance from the critical points: (a) The residual current 
J{p) — J' vs. p — p' for p' values around 0.53685 (the legend 
shows the p' values used) . Due to the finite system, there is 
a small current J' = J(0. 53685), and in order to improve the 
scaling behavior we subtract J' and find that the best scal- 
ing is indeed obtained for p' — 0.53685. (b) The absorbtion 
(=1 — J) for various pc values, cf. the legend. We find that 
the best scaling is obtained for pc = 0.7185. 



